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Abstract 

The problem of structural identification of static plants with 
a vectorial input which contain single-valued nonlinearities, 
is considered. For a problem solution passage in special 
structural space is fulfilled. The method of construction of 
space is described. The criterion of an estimation of degree of 
linearity of plant is introduced. Algorithms of decision- 
making on a class of nonlinearity on the basis of a method of 
secants and construction of area which belong to nonlineari- 
ty are offered. Results of modelling are reduced. 
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Introduction 

Various approaches and methods have been applied 
to identification of nonlinear static systems. In 
[Raybman and Chadeev; Ljung; Graupe] the statistical 
methods based on the correlation and dispersive 
analysis are used. Adaptive methods of identification 
are considered in [Karabutov (2)]. In identification sys- 
tems genetic and neuronet algorithms [Norgaard, 
Ravn, Poulsen and Hansen; Madar, Abonyi and Szei- 
fert] and their combinations with various methods of 
approximation [Righeto, Grassi and Pereira; Sato] 
have been applied. In [Masri, Caffrey, Caughey and 
Smyth; Aguirre, Barroso, Saldanha and Mendes] a pri- 
ori information and the subsequent approximation on 
the set class of functions are presented [Graupe; Espi- 
noza, Suykens and De Moor]. For approximation of 
nonlinear systems in [Billings and Hua-Liang Wei] 
new classes of the wavelet functions minimizing num- 
ber of estimated parameters are employed. Thus the 
choice of a class of functions has not yet been proved. 

Difficulties in structural identification of static objects 
are explained as following factors: 



i) the system exit is the integrated size reflecting influ- 
ence of set input variables; 

ii) absence of the methods, making the allocation of 
necessary interrelations "input-output" for classifica- 
tion of a form of nonlinearity. Approaches to an esti- 
mation of degree of nonlinearity of system are stated 
in [Raybman and Chadeev; Ljung]. The solution to the 
problem of an estimation of degree of nonlinearity is a 
key to a choice of a class of nonlinear static models. It 
is based on application of complex mathematical ap- 
paratus and approaches to structural identification are 
absent; 

iii) attempts of a priori task of nonlinear structure 
[Graupe; Mosteller] on a set of existing inputs in ap- 
plication of the parametrical approach demands the 
problem permission collinearity (multicollinearity) 
[Johnston; Draper and Smith]. The problem is complex, 
bat its decision simple which is defined by conditions 
of practical realization of models. This simple decision 
(exception of dependent variables) is applied widely 
in real control systems. But such approach does not 
allow solving a problem of structural identification. 
Unfortunately, such interrelation is not always unders- 
tood and in existing methods practically is not consi- 
dered; 

iiii) application of parametrical methods on the set 
class of polynoms [Graupe; Ivachnenko and Muller]. 
Efficiency of such approach depends on experience 
and intuition of the researcher. It demands perfor- 
mance of preliminary labour-consuming researches. 
Such approach does not allow defining nonlinearity 
structure in an explicit form. Wide application of pa- 
rametrical methods in problems of structural identifi- 
cation explain their simplicity, which in the given spe- 
cific area fails to give effective results. 

Consideration on specified difficulties, in [Karabutov 
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(4)] the methodology of structural identification of 
nonlinear static systems with a vector input in the 
conditions of uncertainty is offered. Systems contain 
single-digit and multiple-value nonlinearity. In [Kara- 
butov (3)] the method of structural identification of 
nonlinear static systems on a class single-digit nonli- 
nearities (CADN) with a vector input is available. Al- 
gorithms of decision-making on an accessory of a non- 
linear part of system power, logarithmic or exponen- 
tial functions are developed. For logarithmic function 
the special class of transformation of an input is pro- 
vided. On its basis the algorithm of decision-making 
on nonlinearity structure is developed. The decision is 
accepted in special structural space [Karabutov (4)]. In 
this space properties of some nonlinear mapping are 
analysed. The approach to structural identification of 
systems about the CADN and the vector input, offered 
in [Karabutov (4)], makes a multicollinearity problem 
ease [Mosteller]. 

In work on the basis of the results obtained from 
[Karabutov (3, 4)], development of the offered ap- 
proach for static system with several single-digit 
nonlinear functions /(«,) is given. The criterion of an 
estimation of nonlinearity of object on a set of secants 
is provided. For nonlinear object, the decision is 
searched in structural space (SS). The method of con- 
struction of SS is acquired. Structures S v k e (virtual por- 
traits) are put into a SS [Karabutov (1)] which reflect 
properties of a nonlinear part of object. Despite 
nonlinearity S' k e , a decisions on structure /(«,) is ac- 
cepted based on the analysis of linear functions y (se- 
cants S v ke ). It is one of advantages of the offered ap- 
proach to structural identification, which has two 
available methods of decision-making concerning 
structure /(«,) based on application of different crite- 
ria. The case of static object with nonlinearities belong- 
ing to a class of power functions is considered. Me- 
thods and algorithms of an estimation of structure of 
object are acquired provided that, the factor of deter- 
mination of a secant of a virtual portrait should belong 
to some set limited area. The method of decision- 
making on function structure f(u ( ) , proceeding to a 

condition of hit of a secant of a portrait S\ e in some 

sector of structural space is described. The method of 
restoration of nonlinear function in initial space of 
measurements is taken into account. 

Problem Statement 

Consider the object described by the equation 



y n =XU n +B T f{U, «) + £,, (1) 

where U„ e R k , y n e R are input and output an object; 
m, n e U n is the limited irregular function possessing 
property limiting absence of degeneracy; A e Q A c K l 
is the vector of the parameters belonging to limited, 
but a priori to unknown area Q. A ; B e R m ; 
n e J N = [0, N] is discrete time; e R is an distur- 
bance, |£„|<oo, 

/(E/, n) = [/, (« P1 .„ ), f 2 (u p2 J,...f m {u pm ,„)]', (2) 
fm ( u i n ) i s a limited single-digit function, p t e [1, m] . 
For (1), (2) the information set is known 

I =l o (y,U) = {y n ,U n , neJ N =[0,N]} (3) 
and mapping corresponding to it 

r :{#„}->{?„} Vne/^, 
describing an observable information portrait 
[Karabutov (4)]. 

It is necessary on the basis of the analysis (3) and T to 
estimate structure of function f(U,n). 

Further object (1), (2) we will name system as S f . 

Approach to Structural Identification 

The method of an estimation of structure in a nonlin- 
ear part (2) systems S f in the conditions of uncer- 
tainty is based on application of the approach, offered 
in [Karabutov (3)]. We will reduce search of the deci- 
sion to realization of some steps. The first step is to 
estimate the degree of linearity of system and then the 
theorem 1 of [Karabutov (3)] is applied. If the decision 
on nonlinearity of system is true, work in the special 
structure (portrait) S l ke , reflects a state of a nonlinear 

part of system of identification. In work [Karabutov (3)] 
f(U,n) there is the single-digit nonlinear function de- 
pending on one independent variable. Application of 
these results here can appear incorrect. In particular, a 
decision-making can be inefficient a decision-making 
on the basis of one criterion. 

Two approaches are available to an estimation of 
structure of function f(U,ri) . The choice of methods of 
identification of structure appreciably depends on a 
nonlinearity class. More low we will consider a case, 
when f(u j n ) belong to the class of power functions 
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/>,,„) e (F = {«,., e/f:^V«e /„, a,. < oo} . 

The first approach is based on a choice of structure of 
a secant / for a virtual portrait S v k e . Parameter a of 
secant y can be defined so that, to provide a factor of 
determination [Mosteller and Tukey] / of some set 
area. 

We also reduce the second approach to a choice of pa- 
rameter a of function f(u i , n ) e <F . Thus we should as- 
sure belonging of a secant y to structure S v k e to some 

set in special structural space. This set we will name 
sector. 

Estimation of Degree of Linearity Object 

Consider restriction of an observable information por- 
trait r"; c=r J „ V; = \,k . For T"> construct a secant 

y(y,u,) = a 0i +a li u in (4) 

by means of a least-squares method, where a 0j , a lt 
are some real numbers. Define a secant y{y,f{) for 
//=/("/,„) e F> l<k - 
Define set of secants 

r{U, y) = {y(y, «, ), y(y, f x ), i = IJc, I > 1 } , 
set on \ o . 

Consider set on (3) subset of secants for 
S(U,y)aT(U,y), 

which we name a field of secants S s = S(U, y) of sys- 
tem S f . 

Designate <t„ =[fCv„,"i,J Y (y„,u 2 J...y (y n ,u k jf . Con- 
sider the equation 

y n =v T ® n , (5) 

where a vector *P e R k can be defined by means of a 
least-squares method. The vector estimation ¥ exists 
proceeding from above assumptions concerning an 
input U n . 

Structural completeness of system S f in the field of 
structures S s is defined based on the following state- 
ment [Karabutov (3)]. 

Theorem 1. Consider a vector of informative variables 
U„ e R k and set of secants S(U, y) for T B . Then the 
system (1) is linear, if 
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* = 2> ( =1, (6) 
where y/ i is i -th element of a vector *P in (5). 

The theorem 1 gives sufficient conditions of degree of 
linearity (nonlinearity) of system S f on the set field of 

structures S s . If the condition (6) is true, that field 
S s = T(U,y) is full and S s is a linear envelope of an 
output of system (1), (2). If (6) it is not true, value 
N s = 1 - x be termed as an estimation of level of 
nonlinearity. 

Let the condition (6) be exclusive from implementation. 
Hence, we have S s cz T(U, y) . Pass through an evalua- 
tion stage of structure of a nonlinear part of system. 

Virtual Portrait for Decision-making on 
Nonlinearity Structure 

Construct auxiliary set which contains the information 
on function f(U,n) in (1) [Karabutov (4)]. Generate a 

variable s n — I T U n , where / e R k is an unit vector. 

Apply model y s n = a s s n . Define parameter a s e R 

from a condition 

argmin(5> -yj 2 =«;. 

The variable e n = y - y n contains the data about 
nonlinear function f{U,n). Generate set 

l e ={e,U n ,neN}. (7) 

The set l e reflects processes in a nonlinear part of ob- 
ject. Such method of construction \ e excludes domina- 
tion of any input object. 

Now we have reduced a problem to identification of 
structure of function f(U,ri) on set l e . Apply the ap- 
proach offered in [Karabutov (4)], to a choice of argu- 
ments of function f{U,n). Choose such u in , which 

gives the maximum value of factor of determination 
between e and u i . 

Mapping T e :\u in ] —»{<?„} does not allow solving a 
problem of structural identification of system S f 
[Karabutov (4)]. Explain it to that e n e/?has irregular 
features and does not satisfy a condition of Geldera- 
Lipshitsa. Therefore for an estimation of structure of a 
nonlinear part S f construct the mapping, adequately 

describing a state of system of identification on set (7). 
Apply the approach described in [Karabutov (1, 4)]. It 
is based on introduction of the variables, received as a 
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result of processing l e . 
Enter a coefficient of structural properties 

e 

K,u„n =fc s (e,w,.,«) = — • 
»,> 

On J N arrange k eu n is increase. Generate set {^j, 
where k v iq =k s {e ,uj,q) , qGJ' N =[0,N]. As to every 
k v Uq corresponding value e q , receive 

K={e' q , k[ q ,qGj v N ). 

Consider on V k mapping T v ek :{^ 9 }-> { e^'} and 

structure Si e corresponding to it. Si e shows change 

of processes in nonlinear system of identification. 
Properties S v k e are researched in [Karabutov (1, 4)]. 

The virtual structure S[ e has universal characteristic 

for single-digit and multiple-valued nonlinearities. 

The method identification of structure of object (1), (2) 
with one power nonlinearity /(«,) on the basis of the 
analysis S' k e is described in [Karabutov (1, 4)]. In our 
case the offered approach demands updating. In par- 
ticular, except the analysis of properties S\. >e it is nec- 
essary to develop special criteria for decision-making 
on structure f(U,n). 

Algorithms of Decision-making on Function 
Structure f(U,n) 

Application of Method Secants 

We will not consider algorithms of decision-making 
on a belonging /(«,„) to the set class of nonlinearities. 

This problem is researched in [Karabutov (3)]. We 
suppose that f(u ln )e<F,l<m. 

Consider set 

it corresponds to the structure S v k , e , described by 
mapping T v e k . Construct for each structure S v k e secant, 
described by a polynom of degree p 

p 

7 Up { e ^ k i) = fi.p = a oj + Z a U { k i ) ' ■ ( 8 ) 

7=1 

Order of a polynom (8) define so, that 

max J-' ->/>*, (9) 



where r 2 ._ is a factor determination between e and 

« Yi, p 

kj. 

Let 

5^=r^{p'). (10) 

We will describe the approach to structural identifica- 
tion f{u in )&T with the task of border for change of 

factor of determination. Consider mapping 

(«,):{ KMdh {<}- ( n ) 

where («,) is a coefficient of structural properties 
of system identification with an input u"' , a i < oo is 
some number. 

T' ek (a,.) corresponds to structure Sl ^a^m. space 
<P = [k] ,e" ) . We will name <P structural space. 

Construct for a S v k ) secant 

f(e\K q (a,)) = f at = flo/ + a u e\K q («,.) (12) 
and define coefficient of determination r 2 ,_ . 

' To, 

Theorem 2. Let for system S f in space <P : i) the struc- 
ture S v k e and a secant f i p (e v ,kj q j described by the 

equation (8) is constructed; ii) the order of a polynom 
(8) satisfies a condition (9); iii) the mapping (11) the 
description of structure S[ is set; iiii) for the 

K,e{ a i) secant (12) is constructed and the factor de- 
termination is defined r\_ . Then /(«, „) is an element 
of structure of system S f , if the « ; function parameter 
f(u jn ) = u"' is derived from a condition 

r]_ -81 <A,. 7 

where A, > is some magnitude, & 2 is derived from a 
condition (10). 

The proof of the theorem 2 is based on application of 
the theorem 3.6 [Karabutov (4)]. 

Apply the theorem 2 to decision-making on function 
structure f(u jn ), j <m. 

Thus, it is shown that, how on the basis of the analysis 
of virtual portraits S k e (or, ) , S v k e to makes the decision 

on system structure S f . 
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Remark. To check correctness of the received decisions 
the theorem 1 is applied. 

Decision-making on the Basis of Sector Construction 

The decision on F(U,ri) we can accept on the basis of 
construction of set (sector) Sec a (e) , belongs to nonlin- 

earity. The method of construction of sector for static 
systems is offered in [Karabutov (1)] and developed in 
[Karabutov (4)]. We will search for sector on a plane 
(e v ,fc, v ) as area 

Sec a (e) = ( f{e, fc e „ ) , y (e, k eSj a )) , 

limited to the specified straight lines. The area Sec a (e) 
should belongs to a secant y a in space <P . 

Describe a method of identification of a vector F(U,ri) 
on the basis of sector construction [Karabutov (4)]. 

Let the u i condition be r u e > 8 e , where 8 e > some set 
size, is satisfied. Consider mapping r e a {k eu .„}x{e„} 

and its restrictions r*' "' Vz = \,k . For every r]'* a se- 
cant is constructed 

r(e,k eiUi ) = yi =< +<A,„, • (13) 

Corresponding secants y a (12) construct for mappings 

r)*i"> Vu" 1 g<f (j> i) 

Introduce set of secants 

T(K,e) = {T(K e ,e),T(K ea ,e)} V e „el e , 

where 

T(K e ,e) = {y(e,k^),i = lk} , K e =[k eU[ ,...,k eu J , 
K», , a ,^K a K ea e R m m < k T(K ea , e) = {y a j > 1 } 

The bottom boundary of sector Sec a (e) is derived 
from the following statement. 

Theorem 3 [Karabutov (4)]. Any secant 

y(e,k eu )GT(K e ,e) at j = i limits an element 

r(e,k eiUjta .) g T(K ea ,e) from below. 

On the basis of the set analysis T(K e ,e) it is impossible 
to define the upper border of sector Sec„ (e) . 

Therefore the received field of structures are added to 
a secant y(e,k esd ), where s^R, s t = I 7 '(U\u,) . Due 

to the upper border of area of sector Sec a (e) we can 

use a secant y(e,k e s d ) [Karabutov (4)], where d i <oo. 



The statement is true at least for 

Statement. If y (e, k e a a ) c Sec a (e) , which is an 
/(m,) e <F element of structure of system S f . 

The statement proof is based on application of the 
theorem 3.7 of [Karabutov (4)]. 

So, methods of structural identification of a nonlinear 
part considered are based on the analysis of virtual 
portraits of system identification in space <P . The basic 
advantage of virtual portraits S v ke is possibility to ap- 
ply methods analyzed on a class of linear functions. 

Restoration of Initial Nonlinear Function of 
Object 

Consider the problem on restoration of function 
/(«, „) on the basis of the above received results. Let 

the decision on structure be accepted as f(u jn )e'F 
and the set I being irregular character. Hence, func- 
tion change /(«,„) also has irregular character. To 
receive a function form /(«,„) in the specified condi- 
tions on a plane (u t ,f), the following actions are exe- 
cuted. 

In space (P define a secant y a (12) for structure 
$1 e{ a i) ■ Apply model (12) and make the prediction 
e v n of an output of system of identification e v n . Gener- 
ate set 

il(«,)={<, K<,( a i)> i sJ n}- 

On a basis ^(a,.) and (12) an co n estimation is re- 
ceived for w"; 



On a basis co n define an 9 j n input estimation m n . Ar- 
range 9 in and co in on increase and receive a regular 
estimation of function /(«, „) - 

As co i n calculated on the basis of the variable predic- 
tion e n , the operation of preliminary processing of 
values co i n is employed. 

Example 

Consider object (1), (2) with A = [l.5;2.5;5] r , 
B = [0.3;0.3] r , F(U) = [u° A ;uf 6 J . is a stochastic 
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variable with a zero expectation and a final dispersion 
[£„ |< 0.3 . Inputs M[ „ e {/„ are normally distributed sto- 
chastic variables with a final expectation and a disper- 
sion. 

Estimation of degree of nonlinearity of object is de- 
fined by means of the theorem 1. Construct a field of 
secants S s = S(U, y) and calculate a vector fefi 3 in (5) 

and x - Y = [0,339; 0,4; 0,29f, / = 1.03. As de- 
fine 5V X = 0.03 . The object is the nonlinear. 

To structure identification F(U) at first, a method of 
secants is utilized. Function / («,) = «," 4 is also in con- 
sideration. For an estimation of structure function 
/(«)) generates set Y k and construct structure S v k e , 

which is shown on Fig. 1. Receive for a S[ { e secant y p 
(8). Results of modelling have shown that p = 3, and 

coefficient of determination (CD) r 2 „ =0.985 . To a 

e ?p 

parameter choice or, in f(u l ) = u"' apply the theorem 2. 
Suppose A, =0.0012. Structure S' k (a^ is shown on 
Fig. 1. 

For a secant y x of a virtual portrait S\ , (or,) at «, = 0.4 
receive r^, =0.9862 . A parameter choice a, has 
been shown on Fig. 2. a, = 0.4 . 









o 






/Of] \^ 


□ 


o " 



-J H < 1 1 1 1 1 1 1 1 1 1 1 1 1 

-2,0 -1,5 -1,0 -0,5 0,0 0,5 1,0 1,5 

FIG. 1 ESTIMATION / (m [ ) BY MEANS OF SECANTS 



' 0,99 







2 \, 

a, \ 



FIG.2.CHOICE «, 



Compare structures 5^ e , ^ f (a,) (Fig. 1). The correct 
choice of an input gives structure straightening 
S[ e { a i) ■ I* speaks from adequacy of a choice of pa- 
rameter «, . From Fig. 2 the specified estimation has 
been acquired for a { which is equal to 0.42. 

Similarly an estimation for a 3 function /(m 3 ) has been 
defined. a 3 =0.62 (Fig. 3). 



1 ,00 




0,2 0,4 0,6 0,8 



FIG. 3 CHOICE a 




FIG. 4 SECTOR FOR F(U) 

Example of construction of sector Sec a 3 (e) is shown 
on in Fig. 4 in which secants f. (8) are designated 
with p = \ . Function restoration /(«,) is shown in a 
Fig. 5. 



2,0 




0,0-1 , 1 , , , 1 , , , 

1 2 3 4 5 u 



FIG. 5 RESTORED FUNCTION /(«,.) 
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Conclusions 

The problem of structural identification of static ob- 
jects with a vector input and several power nonlineari- 
ties in the conditions of uncertainty is considered. The 
method of an estimation of degree of nonlinearity of 
object on the basis of the analysis of completeness of a 
set of secants of an observable information portrait is 
offered. The virtual portrait (VP) states of nonlinear 
system of identification in special structural space are 
introduced. The method of decision-making on a class 
of nonlinear functions based on the analysis of proper- 
ties VP is developed. The method of identification of a 
nonlinear part of object on based on the construction 
of sector, which belongs to nonlinearity is described. 
Results of modelling have confirmed the working ca- 
pacity of the offered methods and algorithms. 
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